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We calculate the evolution of the geometric entanglement entropy following a 
local quench in the D1D5 conformal field theory, a two-dimensional theory that 
describes a particular bound state of Dl and D5 branes. The quench corresponds to 
a localized insertion of the exactly marginal operator that deforms the field theory 
off of the orbifold (free) point in its moduli space. This deformation ultimately 
leads to thermalization of the system. We find an exact analytic expression for the 
entanglement entropy of any spatial interval as a function of time after the quench 
and analyze its properties. This process is holographically dual to one stage in the 
formation of a stringy black hole. 



Consider an initial, smooth configuration of matter that collapses into a black hole. 
There are longstanding questions about how the information in the initial configuration, 
such as the entanglement between various subsystems, becomes encoded in the result- 
ing black hole. As a quantum theory of gravity, string theory addresses many of these 
questions. While the AdS/CFT correspondence leads immediately to the proposal that 
certain black holes are dual to thermal mixed states of a dual conformal field theory (CFT) 
BU, this says little about the formation process. To go further toward answering such dy- 
namical questions, one needs to study the unitary evolution of a CFT with a gravitational 
dual, undergoing thermalization. In this paper we begin such an investigation in the 
D1D5 CFT, which describes a bound state of Dl and D5 branes, and is well-known as a 
useful system for studying black holes in string theory. 

We study this process using the entanglement entropy, defined below, which measures 
entanglement between subsystems in a quantum system. As opposed to the many studies 
of black hole entropy as entanglement entropy of different parts of the bulk spacetime or 
between different boundary CFTs (0 reviews many of these), we are considering the 
evolution of the entanglement entropy of subsystems of a single CFT, in order to study 
its thermalization. The connection between quantum entanglement and thermodynamics 
has a long history. See [3] for a short review, ^ 13 for relevant early investigations and 
ffoT-tTOj for recent general results. 

As a brief review, begin with a quantum system with Hilbert space H and Hamiltonian 
H. Then factorize, or coarse-grain, the Hilbert space of the full system as H = Ha ®'tIe, 
where the first factor contains all states describing degrees of freedom in a subsystem 
A and the second factor all states for the exterior of A — the "environment" E. We call 
a quantum system, in a state specified by a density matrix p, thermalized if for general 
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small subsystems A the reduced density matrices pa, obtained by tracing over He, are 
approximately (Gibbs/ canonical ensemble) thermal mixed states, e.g., 



e -PH A 



PA « (1.1) 

for some (3 (the inverse temperature), where Ha = Tte H and z 7 ^ = Tr^ e~^ HA . This is an 
admittedly imprecise definition of "thermalized," but is sufficient for our purposes. The 
expression in is appropriate for those cases where the Hamiltonian is the only con- 
served quantity. For systems with more symmetries, including integrable systems with 
infinitely many conserved quantities, one can still define generalized Gibbs ensembles 
and appropriate thermal reduced density matrices, characterized by generalized chemi- 
cal potentials in addition to /3 ©H]. This is relevant for CFTs like the one we consider, 
which undergo a quench but which subsequently evolve as free theories and are charac- 
terized by a set of momentum-dependent temperatures IIT0l[TT| . 

We are concerned here with the case that the full system AE is in a pure state, i.e., a 
closed quantum system. Then the fact that pa is mixed comes entirely from the entan- 
glement of A with E. This entanglement can be measured in a variety of ways, but in 



this paper we consider the Renyi and von Neumann entropies, which are given by (3.3) 
below. We choose these quantities for a number of reasons, including nice analytic prop- 
erties and calculability, to be discussed in the paper. Here we just emphasize that they 
let one track the thermalization of various subsystems as well as deviations from strictly 
thermal behavior (see [12] §8.2 for an introductory discussion of this topic). The entan- 
glement entropies, as a function of subsystem, time, and other parameters, can also yield 
much other information about the system. This is explored in voluminous recent work 
in condensed matter physics (see [13] for several recent reviews). 

One can use these quantities, in principle, to investigate the recently conjectured ther- 
malization time for black holes that saturates a causality bound |fT4f, although we don't 
get that far in this paper. As we discuss below, it is technically difficult to quantitatively 
compare the nonequilibrium dynamics we study here to the system in equilibrium and 
the full thermalization process. However, we can still learn a lot from the results we 
present. 

Motivated by the rapid thermalization observed in heavy ion collisions as well as 
the theoretical questions already mentioned, there are many investigations that use 
AdS/CFT to study strongly coupled CFTs far from equilibrium or undergoing ther- 
malization, e.g., |[T5l - l2T| . Some of these ||22l - l25| also use entanglement entropy via the 
holographic entanglement entropy proposal ||26l l27fl (see also ||28| for a recent study of 
the holographic Renyi entropies). These latter investigations use the dual classical ge- 
ometry and so cannot address some of the most puzzling questions about information in 
black holes, which involve quantum mechanics of the bulk theory in an essential way. 

We study the D1D5 CFT at weak coupling, which has long been used in string theory 
to study black holes H29l |30|. Early studies focused on the extremal, zero- temperature 
configurations, whereas we consider exciting to a state that is far from extremal. Recent 
work has studied Hawking radiation in this system in detail [|3TI - l35| and deformations of 
the CFT away from the orbifold (free) point in the moduli space p6l - l38| (see also p9ll40| 
for leading order calculations away from the orbifold point). Because we work at weak 
coupling, i.e. near the orbifold point, we are not in the regime where supergravity is 
a good approximation. Nonetheless the above work indicates that this regime contains 
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much information about black holes. Additional support for this includes several pre- 
cise matchings between gravitational calculations and calculations from the free theory 
data ||4H - |43l . More evidence for the surprising efficacy of the orbifold CFT in describing 
black holes comes from the very recent paper [|44|. There are also general arguments for 
thermalization of CFTs with gravity duals in the large N limit at any finite value of the 
coupling [|45| corroborated by exact calculations in simplified models ||46l l47| . The re- 
sulting weakly coupled thermal state is dual to a "stringy black hole," in that the string 
length is large compared to the size of the black hole (see [45] for further discussion of 
stringy black holes). The above work indicates that such a state should also tell us about 
traditional black holes. 



In particular, the results illustrated in £ VI have a natural explanation in terms of free 
excitations traveling at the speed of light around the S 1 . This picture is similar to the 
CFT description in [|48l |49| of near-extremal supergravity excitations, which, when the 
decoupling limit is relaxed, can periodically escape the AdS throat to the asymptotic flat 
space. The period and rate of emission were reproduced from the same kind of CFT 
dynamics we observe. Thus our results correctly capture some qualitative aspects of 
the supergravity description. On the other hand, we expect that large energy (far from 
extremal) supergravity excitations can back-react to form black holes. In the CFT, this 
corresponds to thermalization. Since the entanglement entropy that we find does not 
persist, but rather has short-time periodic dynamics, we conclude that, as expected, the 
orbifold CFT does not capture this important process. We hope to address this issue more 
quantitatively in future works. 

A closely related precursor to our work is [50], which also calculates the evolution 
of entanglement entropy in a weakly coupled CFT. The authors proposed this as a way 
to study quantum black hole formation and emphasized that the entanglement entropy 
can be thought of as a coarse-grained thermodynamic entropy. However, their CFT (a 
single fermion) has no clear dual black hole interpretation, although it does illustrate 
some general features of the kind of problem we are considering. 

The general process of thermalization in weakly coupled theories is well studied, but 
the D1D5 system exhibits some novel features. In particular, we consider dynamics aris- 
ing from a localized insertion of a particular marginal deformation of the orbifold CFT 
that acts as a local quench, to be described below. The calculation of the entropy produced 
by this quench is the main result of this paper. This is the basic process by which en- 
tropy is generated. [Note added in proof: After this article was written we learned of 
[51J, which contains calculations of similar quantities and appears consistent with our 
results.] 

As emphasized in [45], the familiar semiclassical dynamics of black holes, including 
the puzzling apparent loss of information, should appear in the limit N — > oo of the dual 
holographic theory. Here that would correspond to the limit of infinitely many Dl and D5 
branes. We do not consider that limit here, rather we consider a finite number of branes 
unitarily evolving toward a thermalized pure state as described above. In this paper we 
just analyze the basic process in that evolution. 

In £jn]we review the D1D5 CFT and the marginal deformations that we study. In £ III 



we set up the calculation of the entanglement entropies in the CFT, including a review 
of the replica trick. Next, in [TV we compute the four-point function that we need to 
calculate the entropies, which we do in *|V We illustrate some of their properties in S VI 
We conclude with a discussion of our results and future directions. 



4 



II. D1D5 REVIEW 

The D1D5 system is realized in IIB string theory compactified or[j]T 4 x S 1 with the 
bound state of Ni Dl-branes wrapping the S 1 and N 5 D5-branes wrapping T 4 x S 1 . We 
take the S 1 to be large compared to the T 4 . The near-horizon limit of the geometry is 
AdS$ xS 3 x T 4 , which is dual to a two-dimensional CFT living on the boundary of AdS^. 

The two-dimensional D1D5 CFT has J\f = (4, 4) supersymmetry with SU(2) L x SU(2) R 
R-symmetry corresponding to the isometry of the S 3 . The two-dimensional base space of 
the CFT is given by the cylindrical boundary of AdS% parametrized by time and the S 1 . 
In addition, we can organize the field content using the S0(4)/ ~ SU(2)i x SU(2) 2 sym- 
metry broken by the compactification on T 4 . One can also fix the total central charge 
c = 6NiN 5 from the algebra of diffeomorphisms that preserve the asymptotic AdS^. 
The CFT has a twenty-dimensional moduli space that corresponds to the near-horizon 
twenty-dimensional moduli space of the IIB supergravity compactification. 

There is a point in moduli space called the "orbifold point," analogous to free SYM 
in AdS5/CFT 4 , where the D1D5 CFT is a sigma model with orbifolded target space, 
(T 4 ) NlNs /Sn 1 n 5 - Just as the dual of free SYM does not have a good geometric description, 
the orbifold CFT is far from points in moduli space that are well described by supergrav- 
ity. We wish to study the effect of certain (4, 4) exactly marginal deformations that move 
the orbifold CFT toward points in moduli space that should have geometric descriptions 
and, in particular, should include black hole physics. Even though we work far from the 
supergravity regime, as discussed in the introduction, we can still capture some black 
hole physics. 

We can think of the orbifold model as N±N 5 copies of a (4, 4) c = 6 CFT. Each copy has 
four real bosons that are vectors of SO(4)/, X\ and their fermionic superpartners. See, 
e-gv 1152 II for details. For computational purposes, we map the real cylinder coordinates, 
t G R and y 6 [0, 2-kR), to dimensionless complex coordinates on the cylinder 

w = t + ia — 1 — y —ir — 9. (2.1) 

R R 

Note that we have also incorporated a Wick rotation in this step. We prefer to perform 
most of the computation in the complex plane by further mapping to coordinates 

z = e w z = e™. (2.2) 

In addition to the local bosonic and fermionic excitations of each copy, the orbifold the- 
ory also has twisted sectors: states which come back to themselves only up to an element 
of the orbifold group S^Ns upon circling the S 1 . The twist operators a n (z) are labeled 
by n-cycles and change the twist sector of the theory. More concretely, consider opera- 
tors O w (z) in the ith copy. In the presence of O"(i2...„)(-2o)/ the operators have boundary 
conditions 

{0^ +1 \z Q + z) i = l,...,n-l 
OW(z + z) i = n . (2.3) 

O^(z + z) t = n + l,...,N 1 N 5 



One may also consider K3 instead of T 4 . 
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Let us emphasize that the twist operators considered here are physical components of the 
orbifold CFT, and should not be confused with twist operators introduced as part of the 
replica trick. 

Following ||36l - l38f , we focus on four of the marginal deformations that involve twist 
operators. These operators are believed to be responsible for thermalization in the D1D5 
CFT. The (4,4) supersymmetric deformations are singlets under SU(2)l x SU(2) r . To 
obtain such a singlet we apply modes of the supercharges to of, where we use 
plus/ minus indices to label elements of SU(2) L doublets and dotted capital Latin indices 
for doublets of SU(2) 2 . In [36J it was shown that we can write the deformation operator(s) 
as 







AbK w o) 



2m A 



w 



dw n-t -\ 
G-(w) 



Wo 



(2.4) 



where the factor of 2 normalizes the operator. The operator a\ is normalized to have unit 
OPE with its conjugate 

' (2.5) 



C2,+ (V)0-2~(^) 



Z' — Z 



This implies that acting on the Ramond vacuum [36J 



OA 



{z)\0- R ) {1) \0- B )W = \0 B ) + O(z). 



(2.6) 



Here |0^) is the spin down Ramond vacuum of the CFT on the doubly wound circle 
produced after the twist. The normalization (2.6) has given us the coefficient unity for 
the first term on the RHS and the 0(z) represent excited states of the CFT on the doubly 
wound circle. 



III. SETUP 



Let us now outline the precise calculation we perform. Since we are interested in 
the dynamics of thermalization or scrambling in the CFT, we quench the system and 
then look at the entanglement entropy of spatial subsystems as a function of time. The 
entanglement entropy of subsystems, as discussed, is a very natural quantity to examine 
when discussing thermalization. Happily, there is already some considerable technology 
for computing the entanglement entropy after both global and local quenches in two- 
dimensional CFTs 

The specific quench we consider is a local insertion of the deformation operator intro- 
duced above. Since this operator is believed to be responsible for thermalization, it seems 
natural to consider the dynamics after its application. Moreover, these results should tie 
in strongly with previous investigations [36-38], which showed that the deformation op- 
erator, in essence, effects a Bogolyubov transformation. For instance, the deformation 
operator, when acting on the vacuum, produces a squeezed state of the form 1136 II 



4(*o) |0^) (1) |0^) (2) =exp 



In 



|0fl> • (3-1) 



In this equation, we only show the left (holomorphic) sector and consider just the twist 
part of the deformation. On the left-hand side we have the two-twist operator acting 
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on the untwisted Ramond vacua, which produces many pairs of bosonic and fermionic 
excitations on the two-twisted Ramond vacuum. Note that the Ramond vacua have an 
SU(2)i x SU(2) R spin structure. The coefficients 7 B and 7 F are functions of z given 
explicitly in II36II . The calculation we propose, then, computes the time dependence of 
the entanglement entropy of this squeezed state; although, we do not use the above form 
explicitly. 

The physical setup is as follows. We apply the deformation operator ( |2.4| > to the Ra- 
mond vacuum at time t = and y = 0. We then look at the time dependence of the 
entanglement entropy of an arbitrary spatial interval. Since we are mostly interested in 
how the entanglement entropy changes due to the quench, we subtract off the entangle- 
ment entropy of the vacuum. We act in the Ramond sector since that is the sector relevant 
for black holes. We will sketch the gravitational picture this corresponds to in the final 
section. 



III.l. Review of the replica trick for computing entanglement entropy 

In the remainder of this section, we set up the calculation of the entanglement entropy 
after the quench. Consider a system S with some subsystem A, and its complement B. 
Recall that the (von Neumann) entanglement entropy of A in S is defined as the von 
Neumann entropy of the reduced density matrix, 

S(A) = - Ty a p A log p A Pa = Tr B ps- (3.2) 

The density matrix p s is the density matrix for the full system S. If, as is true throughout 
this paper, the total system S is in a pure state then p s = (ip\. For our calculation, 
the subsystem A corresponds to degrees of freedom living on some interval of S 1 . This 
definition has a number of nice properties that make it the natural measure of entangle- 
ment including positive definiteness, strong subadditivity, and S(A) = S(B) for a pure 
state. In fact, this is essentially the unique measure of entanglement satisfying the above 
properties ||55|. 

Computing the von Neumann entanglement entropy is computationally difficult be- 
cause of the log, so instead we follow H5S1 - I58I and use the replica trick: we first compute 
the Renyi entropy of order n and then analytically continue to the von Neumann entropy. 
Recall that the Renyi entropy of order n is defined as 

S n {A) = -^—\og{Tx A p\) and S vN (A) = lim S n {A). (3.3) 

1 — 77, 

The Renyi entropies are an interesting measure of entanglement even before taking the 
limit to the von Neumann entanglement entropy. In particular, they serve as a lower 
bound on SVn and vanish on an unentangled state. 

Before showing how to compute the Renyi entropy, we first review how to write the 
density matrix p s as a path integral. From there we can easily compute Tr p\ as a path 
integral with twisted boundary conditions. Let us work in some basis with states that 
we will write as \(p); it is perhaps most natural to think of these as shape states (field 
eigenstates), but any basis works. Then, the y?i-y?2 element of the density matrix at time 
T can be written as 

(<p 2 \p(T)\<p 1 ) = (<p 2 W(T))(il;(T)\ip 1 ) 
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{Me itlT Wi)Me- itiT \^) 



(3.4) 



where we have suggestively switched the order of the two amplitudes for reasons that 
should become clear. The state \ip ) is the state at t = 0, which for us is the state imme- 
diately after the quench. We have the product of two amplitudes, each of which can be 
written as a separate path integral; however, it is more fruitful to think of this as one path 
integral with discontinuous intermediate boundary conditions. More specifically, 



{ip 2 \p{T)\ip x )= I D<f>(t) 
D<f>(t) 



(j>(0)=<p 2 



if° T dtL{<f>(t)) 



e J-t 



-T)=Vo 
e if c dt L{4>{t)) 



Dm 



-T)=i> 



e i/- J dtz,(<Kt)) 



<K0)=¥>i 



(3.5) 



BCs 



where "BCs" in the last line indicates the boundary conditions from the previous line, 
and the contour C starts at t = —T goes to t — and then backwards to t = —T. We 
see, then, that we can think of the density matrix as a path integral which accepts two 
boundary conditions at t — 0. We have translated the ipo boundary condition down to 
t = —T to match with previous calculations Il53"ll54"l . 

While this is formally correct, there are a couple of subtleties to address. First of all, 
we should clarify what we mean by the above path integral, since as written we need 
double-valued fields. It is more precise to parametrize C as 



s e [-T, o] 
-s se(0,T) 



se[-T,T], 



(3.6) 



in which case the action in the path integral becomes 

dtL{<f){t))= I dsL((f>(s))- I dsL T (0(s)), 

J-T Jo 



(3.7) 



c 



and is single-valued on s. We put the superscipt T on L in the second term to indicate 
that it is the time-reversed Lagrangian. The second issue we need to address is Wick- 
rotating the path integral. We usually Wick-rotate the path integral to imaginary time to 
make the oscillatory term iS into a convergent — Se', however, we now have a minus sign 
between the two terms in (13. 7), which means that we should Wick-rotate the two terms 



oppositely, see Fig. [T] When we Wick-rotate the second part in the opposite direction, 
we get rid of the minus sign and the time-reversal: we get a smoothly defined Euclidean 
path integral 



Z(to,t/;Vo;¥>i,¥>2) = 
BCs: 0(r o ) = 0(r^ 



D<P(t) 



exp 



BCs 

0(0" 



¥>2, 



drL E (4>(T)) 



(3.8) 



We can compute this path integral for r < < tj with real r and r/, and finally analyt- 
ically continue to the desired matrix element via 



(ip 2 \ P \(fi) = lim N e Z(-iT - e, -%T + e; ip ; y?i, (p 2 ) 

e->0+ 



(3.9) 
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FIG. 1. The contour along the real axis of the complex r-plane for the Euclidean path integral. 
The analytic continuation back to Lorentzian time is shown in light gray, which shows how the e 
regularization arises. 

where we put in e to "remember" which direction we Wick-rotated the two terms. The 
factor of N e is a normalization constant that ensures Tr p — 1, 

±r = (Vol e" (2£) ^ |Vo> = J D <t> Z{—iT - e, -iT + e; (p, (p). (3.10) 

The limit as e — > is both delicate and crucial to getting the right physics, since 0(r) has 
a branch cut along the negative imaginary axis. Later, it should become clear that e plays 
the role of a UV cutoff. 

Before continuing, let us remark that the above should be reminiscent of the Schwinger- 
Keldysh, or closed time path, formalism with temperature T = 1/ (2e) (see, e.g., ||59| for 
a review of this formalism). Indeed, if one integrates over ip , then it is exactly the 
Schwinger-Keldysh formalism, with some insertions at t = 0. Also note that if one iden- 
tifies if i = p 2 = V 9 an d integrates over </?, then one computes Tr p, which is unity for a 
pure state. 

We now have all of the tools to understand how to compute the Renyi entanglement 
entropy as a function of time after the quench. First note that it should now be clear how 
to compute the reduced density matrix ( |3.2| >: 

(a 2 \ Pa = (a 2 | Tr B p |ai> 

= / DbN e Z(-iT -e, -iT + e;ip ;<Pi = {oi,&},y?2 = {a 2 ,b}) 
Jb 

= N e Z A [ - iT - e, -iT + e; ijj ; a 1} a 2 ). (3.11) 

Here we indicate a field taking values a on A and b on B by {a, b}. We can compute 
this quantity from the same path integral in ( 3.8| > with altered boundary conditions at 



t = 0. In the region B, we now demand that <p be continuous at t = 0. We started with a 
full cut, which we sew together in region B. The boundary conditions on the remainder 
determine the matrix element computed. This manifold is pictured in Figure |2j 

Now, to compute Tr^ p\ we start by inserting J Da \a) (a\ in between each pa and then 
perform the trace in the |a) basis. This becomes n distinct copies of the above path integral 
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FIG. 2. The reduced density matrix as a path integral. Note that the flat piece on top is there 
for illustrative purposes only. The boundary conditions on the bottom two edges are both \ipo), 
whereas the boundaries in region A are "inputs" that determine the matrix element of the reduced 
density matrix. 

with appropriate integrals over the a^. 

Tr Pa = J Da Q J Dai - J Da n -i (a \ p A \a n -x) • • • {a 2 \ Pa K) (oi| Pa K) 

= J Da --- J Da n - X (N e ) n Z a (t , r f ; On-i, do) • • • Z A {r , r f ; ip ; a , ai), (3.12) 

where r and Tf get analytically continued as described. One can then put all of the pieces 
into a single path integral over n replicas, with n-twisted boundary conditions in region 
A connecting the replicas and singly-twisted boundary conditions outside of A. 



III.2. Entanglement entropy in the D1D5 CFT 

Let us apply the above general discussion to the matter at hand. We need to compute 
the twisted path integral described above. We can rewrite the path integral as a correlator 
of local twist operators that induce the appropriate monodromy, and then compute the 
correlator using techniques in 1160 H . Let us note that there is an extra layer of obfuscation 
beyond computations in other CFTs since our quench involves a distinct twist operator 
that is part of the physical spectrum of the CFT. 

We prepare the state |^ ) by starting with the vacuum at r = — oo, evolving forward to 
t where we insert our quench O(w ) from ( |2.4| >. To compute the Renyi entropy, we need 
n replicas of |V>o). The trace in ( 3.12| > is then proportional to the four-point function 



w n {r ,r f ;e 1 ,e 2 ) = {[ot( W/ )]" ( r ti (t ( ; 2 ) ( T n K)[OK)]"), (3.13) 

where 

w = t Wf — Tf W\ = 161 w 2 = i0 2 . (3.14) 

Since each of the O's has a 2-twist it becomes necessary to clarify the branching structure 
of the correlator. Let us specify the indices of the twist fields involved in the correlator: 

[a 2 ] n = CT(i2)CT(34) ■ ■ ■ 0"(2n-l,2n) &n = C(135...2n-1) • (3.15) 
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The indices are labels for the 2n sheets involved in the correlator, and we use the par- 
enthetical notation for single cycles of S n . Note that the bare twist a n introduced by the 
replica trick twists one copy from each pair of twisted replicas. This fixes the topology of 
the correlator. 

We can now write the Renyi entanglement entropy as 

W n (-iT-e, -jT + qgiA) \ 
[W x {-iT-e,-iT + e)Y )' 1 ' ' 

Note that a± is the identity operator and so there is no need to specify the Q\ and 9 2 for W\. 
Also note that any normalization issues from defining \ip ) in terms of the local operator 
O(w ) cancel out between the numerator and denominator. 



S n (T, 0i, 9 2 



1 — n 



lot 



IV. THE FOUR-POINT FUNCTION 



The four-point function in Equation (3.13) factorizes into a four-point function of bare 
twist operators that we compute by mapping to a covering space and a correlator of 
insertions in the covering space. 

We first compute the correlator of the bare twists, and then treat the non-twist su- 
percharge insertions that appear in the covering space. We map the correlator in Equa- 
tion ( |3.13| > to the plane via the exponential map \22) . We will then treat the associated 
Jacobian factors in |[V4 



IV.l. The twist correlator 

Let us begin, then, with just the twist part of the correlator 

([a 2 (z f )} n a n (z 3 )a n (z 2 )[a 2 (z )} n ) . (4.1) 

Note that this part of the correlator applies to any CFT with two copies that are suddenly 
joined by o 2 . Thus when discussing the bare twist results, we keep c the central charge of 
a single copy. The SL(2, C) symmetry determines the form of the 4-pt function up to an 
arbitrary function of the cross-ratio. Therefore, we can compute the 4-pt function 

JF n (u) = ([a 2 (oo)]V n ( M K(l)[a 2 (0)n , (4.2) 



and then find the 4-pt function of interest in Equation ( 4.1| >. 



To compute the four-point function we need to find a map to the covering space, and 
then compute the Liouville action associated with the map [60J. Fortunately, Appendix D 
of II52II gives an explicit formula for spherical genus correlation functions of S^-twist 
operators as a function of the coefficients of the map. Once we find the map, we can 
make use of the formula to avoid computing the Liouville action directly. 

Let us now list the properties the map z = z{t) from the 2 -plane to the t-plane must 



have, as determined by the index structure shown in Equation (3.15). First, one can show 



from the Riemann-Hurwitz formula that the covering space must have spherical genus: 



It, 



s + 1 = -[1 • n + (n - 1) + {n - 1) + 1 • n\ - (2n) + 1 = 0, (4.3) 
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where is the ramification of the ith point with nontrivial monodromy and s is the total 
number of sheets (or indices) involved. Second, the map must have monodromy at z = 0, 
1, u, and oo appropriate for their respective twist operators. For example, the point 2 = 
must have n images in the covering space, each with monodromy 2. Third, generic points 
in the 2-plane should have 2n distinct images in the t-plane. Thus, we are looking for a 
meromorphic function z(t) with the following local properties]^] 



z 



a* {t - t*Y z ps 0, t ps t* j = 0, . . . ,n - 1 



z-?i~ a M (t — 1)' 
2 ps 6 i 2 

z 



z ps 
2 ps 



M« 

M, t PS 1 



2 ->• oo,t ->■ 



OO 



(4.4) 



2 — ► OO, t PS t 



J 



n 



where the a's and b's are coefficients that are determined from the map and the t* and t°° 
are the various preimages of z = and z = oo, respectively. We have fixed the points 0, 1, 
and oo in the t-plane. The remaining points in the t-plane must be determined from the 
map. 

A rational function that satisfies the above properties is given by 



At n - (t - l) n 
t n - (t - l) n 



where 
Note that 

and we can write 



u — A . 



5 



1 — Ane 1 ™ 



'.7 



1-e 



; 2irj 5 



71-1 



dz 2U{A-If _ yr _ 



and 



r - (t - 1) 



n— 1 



n-l 



3 ' 



At r 



(t-ir=(A-i)H(t-t*) 

j=0 



(4.5) 

(4.6) 
(4.7) 

(4.8) 
(4.9) 



This allows us to write the map in a form more conducive to finding the a's and b's, 



{a -if n?=o(*-q) 2 
n2 u]zi(t-t^) 2 ' 



(4.10) 



2 Let us note that in an unfortunate notational choice the stars do not indicate complex conjugation. They 
are merely labels. 
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The coefficients can be written as 



a„ 



2(-l) n+1 (A- 1) 
2A(A-1) 

n(A-l) 2 

mr - (n - m 



n-l 



ir 1 n 



bo 



(A -If 



rr 



(A -If 



n- 



TP (t c 

We can plug these into the formula from [|52| to find F^u) 



(4.11a) 
(4.11b) 

(4.11c) 

(4.11d) 
(4.11e) 



M 



Fn{u) = n 



Pi 



;(Pi+l) 



v i=l 




9i+l 



12 g- 



where 



P.7 



Pi = P„ 



n 



l & o| g g| ; 

(4.12) 

(4.13) 
: 6 for the 



and c is the central charge of a single copy of the CFT. (Remember that c 
D1D5 CFT.) All of the various coefficients and products can be checked by appropriate 
Laurent expansions and by judicious use of the identities ( 4.9| >j^| The final result takes a 
very simple form 



F n {A) = (2-f"n-«(" +1 )) (2i^) \AA(A- lff^ 



2n 



A 



2n-2 



n 



{A -I) 



2n-4 



(A - lf n 


8 


(A -If 


n 2n+A 




n 2 



2^ 



\AA(A-1) 



An- 1 ) 



(4.14) 



where recall that A = y/u. As usual, this statement is ambiguous since the square root has 
two branches. For instance, as one moves u to z — 1, A approaches either +1 or —1; the 
amplitude is zero in one case and not in the other. The location of the branch cut plays a 
crucial role in getting the correct answer, so let us take a moment to discuss its physical 
origin. 

There are two sets of branch cuts associated with the correlator in Equation (4.2 ). There 
are the (12) (34) • • • branch cuts which extend from the origin to infinity, and there is the 
(135 ■ ■ ■ ) branch cut that connects the two cr n 's. From this structure, we see that each of 
the two a n 's can be on one of two branches. Physically, then, when u approaches 1 the 
two <r n 's could be on the same branch and colliding, or they could be on two separate 
branches and widely separated. This explains the physical origin of the square root. 



3 One needs to plug in with specific values of t and take derivatives to get the desired results. 
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IV.2. The non-twist insertions 



We now include the effect of the non-twist operators. Each application of G_i brings 

in a Jacobian factor when mapped to the covering space. The chiral primary inserts 
an appropriately normalized spin field S ++ in the covering space [61J; all other Jacobian 
factors associated with the chiral primary twist operator are taken care of by the Liouville 
action and the operator normalization. 

Thus for each operator inserted at the origin, we write (this is the holomorphic part) 



dz 

2^' 



z)<r+(0) 



V2 



dt fdz\ 
\dt) 
dt 1 



'dx 



Jt * 2m y/t-t* 
-dX AA (t*)S A (t*). 



AA 



(t*)S A (t* 



(4.15) 



The factor of a* in the square brackets with S + is a local normalization that comes with 



the definition of crj Q6TJI- We have suppressed the j index on a* and t* in the above. Let 



us define 

A (t)=dX AA (t)S A (t). (4.16) 

For the operators inserted at infinity we can write a similar expression; however, it 
is helpful to make what we mean by "infinity" more precise by making the complex 
plane into a sphere. Following ||60l|, we cut both the z- and t-planes into large discs with 
all operators in the finite plane enclosed, and then glue equal-sized discs on top. Then 
infinity of the original plane becomes a point centered on this "second" disc. The radius 
of the 2 -plane discs is 1/8 and the radius of the t-plane discs is 1/8'. 

The insertions at infinity become 

(G t ) A _ | a 2 , + (oo)) 8i bl (OY(n, (4.17) 

where again we have suppressed the index j on b and t°°. The exception to the above is 
the j = insertion at = oo. This insertion, gets an additional factor of 8': 



(G*)i _x<7 2)+ (oo)) *i bl {OY{t = 0). 



(4.18) 



The above equations were written for just the left part of the operators, but we must 
also include the right part. Thus, the contribution from the non-twist insertions may be 
written as 



5^ 



: (OY B (ttKOY B (t?) ■ ■ ■ (0V*(Ci) : = O AB (t* ) ■ ■ ■ AB (t* n 



(no sum) , 



(4.19) 
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where the insertion at should be thought of as at t = 1/8' (really at t = 0); the end 
result being that the S's cancel out in the limit as 5' — > 0. We can evaluate the products of 
the a*s and bs as before: 



5 5 2 n 



n 



\A — lj 5 ^- 1 ) 

5(n+l)M|5^ 



no sum . 



(4.20) 



The bracketed expression is the Jacobian factor from the map and the factor of <5 5n / 2 cor- 
responds to the fact that we have n operators with weight 5/8 (subtracting off the weight 
of the bare twist) in the left and right sectors inserted at infinity in the z-plane. Note that 
the correlator of bare twist operators has a factor of 5 3n / 2 that has been dropped in ( 4.121 ), 
since we really want the regularized correlator. 

Since we are dealing with free fields the above correlator is easily evaluated in terms 
of Wick contractions. For instance, for n = 2 we have 



-AOO J.* 

ri t o 



*|5 



Uoo j.* 1 5 



1 1 


-5 

+ 


1 1 


-5 


2 1-VA 




2 1 + y/A 





m-Va\ 10 + \i + Va\ 10 



a| 5 



(4.21) 



Thus putting all of the contributions together for n = 2 we find 



T 2 = UA(A-lf 



23 . . 

2 \A\ 



13 . 

4 L4 



\A-1\ 
2 15 \A\l 

1 



\i-Va 


10 + | 


i + Va 


10 




1-A 


5 



i|-I (\i-Va\ 10 + \i + Va 



10 



(4.22) 



We use this to find the Renyi entropy of order 2. 

Unfortunately, we could not find a closed-form expression for the general covering 
space amplitude suitable for analytic continuation. It may, however, be written as the 
sum over all total Wick contractions: 



(: (OY B (oo)(OY B (t?) ■ ■ ■ (Ot)^(d) : : O Ab {Q ■ ■ ■ O^(Ci) = ) 

n-l 

SGSn 3=0 



where t j>k = tf -t* k . 
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IV.3. The four-point function of interest 



We can use (4.14) to find the 4-point function that we actually want. The general form 
of the 4-point function is dictated by SL(2, C) symmetry to be (cf. II62II ) 

a, = (M^M*s)Mz*)M*i)) = f(v, fj) n 4 "'4 "' v = ( 4 - 24 ) 

Z13Z24 

where z^ — Z{ — Zj and f(r), fj) is a function that is completely undetermined by SL(2, C) 
symmetry We can compute the 4-point function with points 0, 1, u, and 00, and determine 
/ and therefore the general 4-point function. 



IV.3.1. The four-point function of bare twists 



We separately compute the entanglement entropy of the bare twist operator a 2 and of 
the full deformation operator. 

In our case the operators are left-right symmetric and therefore hi = h if and 



Aa 



m I! 



-u 1 



|2(§-/ii-/ij) 



i<3 



We know that the conformal scaling dimensions for the bare twists are 



n- 



24 



1 



Above, we computed the correlator with 

Z\ = z 4 = 00 

For this case, we have 

V = 

and thus the 4-point function takes the form 



ho = h* 



z-i 



c 

24 



n 



n 



z 3 = u. 



u 



00 



-4/14 



f{ri)\l-u\ 2 ^- 2h ^\u\ 2 ^~ hl ~ h2) 



(4.25) 

(4.26) 

(4.27) 
(4.28) 

(4.29) 



where we regulated the factor of "00" by putting the CFT on a disc (see above). We are 
interested in the finite part. 

We can then write the general 4-pt function in terms of the one we computed as (we 
use the fact that hi = h 4 and h 2 = h 3 ) 

-4/12 



Aa 



JFJu) 



V 



V 



-Z23 



l*b/l 



-4fti 



U 



1 

v' 



Plugging in with the above weights and with T from Equation (4.14 1 



(a 2 (zf)a n (z 2 )(T n (z 3 )a2(zo)) = \z } - z \ 



(1 + VvY 



where recall 



(zq - z 2 )(z 3 - z f ) 



16rj(z 2 - z 3 y 
(z - z f )(z 2 - z 3 ) 



(z - z 3 )(z 2 - z f ) (zq - z 3 )(z 2 - z f ) ' 

Note that the amplitude is invariant under interchange of z 2 and z 3 or z and Zf. 



(4.30) 

(4.31) 
(4.32) 
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IV.3.2. The four-point function with the full deformation operator for n = 2 



For the full deformation operator the weights are given by 



hi = /i 4 = n h 2 = h 3 = ■ 
Plugging in as before, this gives 

\d^z f )} n a n (z 3 )a n (z 2 )[d(z )} n ) = F n {u) 



n 



n 



%23 



-(n-i) 



- i» 



where as before we should replace u with 1/rj. 



For the case n = 2, we can plug in with T 2 from Equation (4.22 1: 



[d\z f )} 2 a 2 (z 3 )a 2 (z 2 )[d(z )} 2 ) = 2~" \z 23 \^\z 0f \ 



(4.33) 



(4.34) 



1 - 7?i| 10 + |1 + Tj- 



110 



(4.35) 

We have carefully written the above expression so that the r\ \-t 1/rj symmetry is manifest. 
This symmetry comes from the exchange symmetry z 2 H Z3 or z f> z/. 



IV.4. From the cylinder to the plane 

The physics of the D1D5 system originates on the (Lorentzian) cylinder, so we should 
be careful to put in Jacobian factors that arise in using the exponential map from the 
cylinder to the plane. Normalized states are normalized states, so we do not need to 
worry about Jacobian factors for the O, and 0\. We do, however, need to worry about 
Jacobian factors from the replica twists. 

We started on the plane with dimensionful coordinates t G R and y £ [0, 2nR) and 



then introduced Euclidean dimensionless coordinates r, 9 via (2.1 ) that may be written as 



Dls wrap. Finally, we map to the complex plane using the map in (2.2 1: 



a complex coordinate w = r + i9. Note that R is the radius of the large S 1 cycle that the 

(4.36) 



z = e w = e T+ie -> e - l Ti +i R 



z = e* = e T - ie -> e-^R, 

where the arrows show how we should analytically continue back to real time at the end 
of the calculation. 

This is a convenient point in the discussion to put in the appropriate normalization so 
that 

N {0 ] 0) — 1 =^ N — Izo/T 1 - (4.37) 

This is in the 2-plane, but as mentioned above, we can put in the normalization on the 
cylinder or on the 2-plane. One can check that the normalization ensures that we insert 
a normalized state, which means that W\ = 1. Note that this factor cancels out the in- 



verse factor in Equation (4.30 1. Starting from the cylinder, we can compute the four-point 
function as 



W n = \z of \^ (Ola n (w)a n (w + iA9)O t 
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l4fei 



4A n / 

e( n -^) / t 



ze 



iA8 



)0, 



0]a n (z)a n (ze^)0 



(4.38) 



This normalization ensures that W\ = 1 and this will let us easily compute the entan- 
glement entropies. Note that while the Jacobian factor is unity for n = 1, it still gives a 
nontrivial contribution to the von Neumann entropy defined as the limit as n — > 1. 



V. THE ENTROPY 



We now have all of the pieces to discuss the entropy. Let us first treat the Renyi entan- 
glement entropy of the bare twist operator: 



log 



W„ 



l-n"° W{ 

c n + 1 

log 

24 n 6 



R 



w 




z 3 


1 


16r] 


\%2 


2 


Z3 


2 


(i + Vv) A 



(5.1) 



Since we have carefully normalized the W n so that W\ = 1, the von Neumann entropy is 

c 



S vN = lim S n 

mi 



12 



log 



R 



w 




^3 


4 


I677 


w 


2 


Z3 


2 


(1 + ^ 



(5.2) 



The limit is trivial only because we already canceled out the 1 — n in the denominator. We 



illustrate various properties of this formula in { VI 

Let us also write down the n = 2 Renyi entropy for the full deformation operator. This 
is given by 



So 



log W 2 = - log 



„_M 3 I^N* 

2 2 _R 2 3' 

1*23 1 5 



1 + v 7 ^ 



1 - 77^ 1 10 + |1 + r/2| 10 

1 n 



(5.3) 



Recall that the Renyi entropy gives a lower bound for the von Neumann entropy and that 
it vanishes if and only if the reduced density matrix is that of a pure state. 



V.l. Entanglement of the vacuum 



It is useful to subtract off the contribution to the entanglement entropy from the vac- 
uum, which we compute here. We want to compute this as a function of the physical 
cylinder coordinates, so we include Jacobian factors for the map z — >■ y, while keeping 
the convenient variables z = e %y l R . From the above, we have 



QV 



1 



1 — n 
1 

1 — n 



log 
log 



Kn-i) 



1 z 2 
Rz 2 - z 3 



a n (z)a n (ze i0 )) 
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c n + 1 
6 n 



\ogR 



Z2 ~ Z 3 
Z2 



(5.4) 



The limit as n — > 1 yields the von Neumann entropy 



S™ c = -\ogR 



Z2 ~ Z Z 



Z-2 



(5.5) 



Since z 3 = e z 2 we have 



z 3 ~ z 2 



Z-2 



Z3 
Z-2 



21 sin 



21' 



(5.6) 



and so 



cvac _ ^ i 
D vN — g 10 S 



~ L ■ i 




— sin 




7T 





(5.7) 



and L = 2tiR to 



where we have written the above expression in terms of Z = \y 2 — y\\ 
show agreement with [57, 58 j for a CFT on a cylinder. 

Actually, in the above we have suppressed an ultraviolet cutoff a in terms of which 
we measure the physical lengths I and L. So the correct expression is (5.7 1 with / and L 
replaced by l/a and L/a. There is a logarithmic divergence of as a -» as noted in 
early investigations |57ll64| . One can either treat the expression as an asymptotic result 
for a small but finite a, as would be appropriate for studying lattice theories (e.g., in 
condensed matter theory), or find "renormalized" entropy differences between various 
states in the theory. In some cases these differences are finite in the a — > limit [57J. 

Here we follow the latter approach (although, as it turns out, our entropy differ- 
ences are still divergent), since we are interested in the extra entanglement added by 
the quenching process: the entanglement entropy increase, AS. After subtracting off the 
entanglement of the vacuum we get for the insertions of the bare twist operators 



ASV N = — log 



16r] 



(5.8) 



and similarly for the full deformation operators: 



ASo = loe 



23 

2^ 



1 + 



|1 - 773 1 10 + |1 + 773 1 10 



(5.9) 



VI. PROPERTIES 

In this section we include plots of the entropies ASVn and AS 2 for a variety of space 
and time intervals and examine some of their properties. We discuss the details of the 
UV cutoff and the choices of branch cuts needed to compute with the entropy formulas 
we have given. We start with the entanglement entropy introduced by the bare twist 
operator, then consider the full deformation operator. 
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VI.l. The entanglement from a bare twist 

Our formula ( |5.8| > can now be computed for any time and choice of interval [9\, 9 2 ). 



We recall that the quench is located at 9 = and t = 0. To actually compute from ( |5.8[ ) 
we must choose c, e (as desc ribed in yil.ll >, and a prescription for the square roots (as 
mentioned at the end of £IV.l I. In this subsection, we choose c = lJ^The entropy increase 



AS" diverges as e — > for some values of t, see (6.9 1 below. We can understand this by 
thinking of e as a UV suppression factor or regulator, consistent with its appearance in 
the evolution operator e~ lH ^~ ie ' . it is perhaps not surprising that the local quench has a 
UV divergence since it is localized at a point, a phenomenon observed in Il65l . In a real 
physical process there should of course be a finite amount of energy, which we can treat 
as a finite e in our calculation. We would fix the appropriate value of e in a full treatment 
of the physical evolution of the D1D5 CFT, but here we leave it as an unfixed but small 
parameter. With finite e we have 

_ sin (§(| + 9 1 - ie)) sin (|(| + 9 2 + ie)) 
sin (|(| - #x - ie)) sin (§(£ - # 2 + ie)) ' 



sin 



(!(!-0i + ie))sin(f(!-0 2 -ie)) 



where here 0j = yi/R and e should really be e/i? so that all parameters are dimensionless. 
For the remainder of our discussion, we set R — 1. In the above expressions and in what 
follows, 9 parametrizes the double circle by running from 9 = to 9 = An. This explains 
the factors of 1/2 inside the trigonometric functions. 

Using sundry trigonometric identities, we can rewrite 77 and f] in a slightly more useful 
form 

cot *±&- cot ^ , x 

77 = e 2 ^ tan cp = e 2 —— (6.2) 

' ^ l + e 2 cot^cot^' V ^ 

where we have introduced e = tanh(e/2). For our purposes, it suffices to drop the e 2 term 
in the denominator and write 

tan Lp w e [cot _ co t « ± => cp w nvr n G Z. (6.3) 

From the above, we conclude that for small e tp is close to some multiple of it, and there- 
fore r\ is essentially unity. This fact could have been read off from ( 6.1| > directly; the key 



realization from the above is that the sign of the depends on time which implies <p has 
some nontrivial time dependence. If we plot cp as function of t for reasonably small e 
and take care with the signs, we get Figure [3] Note that tan<^> never vanishes in ( 6.2[ >, 



and therefore — n < cp < 0. We have chosen these particular branches to be consistent 



with (6.41. 



All of the time dependence comes from the interaction of this phase with the square 
root. We have carefully chosen ip and Cp in Figure [3] so that the correct branches are 



^r\ = e lLp a/77 = 772 = S = - e *§, (6.4) 

4 The appropriate choice for the D1D5 CFT would be c = 6; however, the c-dependence is just an overall 
constant here, and since the bare twist operator is not the deformation operator, we may as well just set 

c = 1. 
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FIG. 3. Plot of ip (blue, dashed) and (p (red, dot-dashed) versus time for e = 10~ 2 , where rf = 
exp 2ip and fj = exp 2iip. Note that ip and ip obey the strict inequality — ir <p<0<ip<ir for 
all time. In the limit as e goes to zero, the function converges pointwise to a piecewise function 
saturating the inequalities. Indeed, for e = 10~ 4 the plot is (visually) indistinguishable from the 
corresponding piecewise function. 

where the fourth roots arise when considering the Renyi entropy of the full deformation 
operator. Then, we can write the entropy in terms of if and (p 

c c 
ASW = - - log ( cos \ip) - - log ( cos \tp) (6.5) 



Note that in defining cp and ip from ( 6.2| > there is an ambiguity associated tan; similarly, in 



writing the square root of rj in terms of <p and (p there are two branches one could choose. 
However, the above choices are fixed for us by causality: the entanglement entropy of 
our interval cannot change from the vacuum value until a signal traveling at the speed 
of light from the quench could reach the interval. The origin of the branch cut is also 



discussed at the end of £IV.l 



Since for small e, ip and (p spend the vast majority of time near or ±7r, let us examine 
the limiting behavior of the entropy away from the transitions. This is where the approx- 
imations in Equation ( |6.3| > are good. Let us define x to be the difference of cotangents, 



x = cot ^ - cot (6.6) 



2 2 

Then, for small ex we can write 

\{lx) 2 + 0((ex) 4 ) x<0 




(ex) 2 + 0((ex) 4 ) x>0' ^ 

and thus the left-moving contribution to the entanglement entropy is given by 

c i I v\ f|(ex) 2 + 0((ex) 4 ) x<0 
Sl = -log ( cos^) w < 8 , /'.o^ ■ (6.8) 

6V 2> \-\og\ix + 0({ixf) x>0 V ' 

We see that the entanglement entropy (away from transition regions) either vanishes like 
e 2 or diverges like — log e. The peak value occurs at t — — ^ 2 + 2nn for integer n. We 
can estimate the peak value as 

S p L eak ps - log e + log (tan ^) + 0(e 2 ). (6.9) 



21 



The unbounded growth of S^ e3k as e — > indicates that the state after the local quench 
has entangled elements of arbitrarily high energy, a UV effect already discussed. 

The e dependence of S^ eak indicates only partial thermalization, as we now discuss. We 
can think of the UV regulator e as introducing a temperature l/2e, in that it corresponds 
to introducing the operator e ~ 2eH in our density matrix given in (3.4). We can compare the 

e dependence of S\ with the entanglement entropy of the interval in the equilibrium 
mixed state at temperature l/2e. The limit e — > corresponds to high temperatures and 
in this limit the entropy should be dominated by the extensive quantity 

AO 



S max 



(6.10) 



This gives the asymptotic behavior of a general CFT of length AO at temperature 1/e, 
which is a regularized expression for the maximum entropy on that interval. Thus, we 
see that while the entanglement entropy we produce diverges, it is parametrically less 
than the maximum possible entropy for the subsystem. In other words, even at peak en- 
tanglement, the system is far from being "Page-scrambled" as defined in lHHl , referring 
to ||66|. This is expected since being Page-scrambled would require the reduced density 
matrix after the quench to evolve to the maximal-entropy (infinite-temperature) thermal 
density matrix, proportional to the identity operator on Ha, but we do not expect a lo- 
cal quench to lead to a thermal reduced density matrix (at any temperature). First, this 
system has decoupled momentum sectors which can thus be independently thermalized 
with various momentum-dependent temperatures [fTT||. Second, the local quench pro- 
duces coherent sets of noninteracting particles traveling from the quench point, so there 
is no mechanism to scramble their momenta. So thermalization must involve more than 
the process we study here, as we discuss further in the Conclusion. 

Now we can examine some specific calculations of ASVn, choosing c = 1 and e = 
10~ 4 . The branch cuts are chosen as discussed above. The first obvious feature from the 
formulas is that the entropy is 2n -periodic, which follows from the 27r-periodicity of the 
quenching process; the point where the two circles join becomes two antipodal points on 
the length AirR circle. 



In Fig. 4a we choose the interval [n/2, 3tt/A] and see positive entropy at those times 
when the null world line from the quench point intersects the interval. We have sepa- 
rated the contribution from the left-moving and right-moving sectors for illustrative pur- 
poses. With the interval to the right of the quench point, the entropy can be qualitatively 
understood in terms of particles emitted from the quench point and traveling with unit 
velocity, a picture first described in |53|. The positive entropy comes from the presence 
of entangled pairs of (left- or right-moving) particles, one member of which is inside the 
interval and the other outside, and so is traced over. A space-time picture of this is given 
in Figure [5] and we discuss it further in the Conclusion. AS 2 behaves the same way, as 



shown in Figure 4b In fact, plots of AS 2 take almost exactly the same shapes and share 
all the qualitative properties of those of AS vN , so we only include plots of the latter in the 
following. 

We can translate the interval toward and away from the quench point, as seen in Fig- 
ure [6a| which is seen to have almost no effect on the entropy curve. However, if the 



quench takes place inside the interval there is a very noticeable effect shown in Figure 6b 
due to overlapping contributions from the left and right-moving sectors. We can examine 
intervals of different sizes, as in Figure [7| and see a clear dependence on the size of the 
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FIG. 4. AS, 



(a) 



1 1 2 JT , i , 3 7T , . 



(b) 



4a 



and 4b for the interval [vr/2, 37r/4]. For ASVn the right-moving contribution 
is shown in red, the left-moving in black. For AS2 we cannot separate the left and right-moving 
contributions. 




FIG. 5. Here we show the two circles being joined at 9 = and 6 = 2ir. Time is increasing 
up on the figure. The quench occurs at t = 0. Light cones are emanating from the two quench 
sites. The right-moving excitations travel along the green curves, and the left-moving excitations 
along the blue curves. The grey vertical strip represents the time evolution of the interval [0i, 62], 
and the red parts of the strip show when we expect nonvanishing entanglement from the particle 
interpretation. 



interval. The dependence of the peak value on the size can be read from the second term 
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(a) 




(b) 



FIG. 6. 

xn/20 
0,1,2, 



ASVn for various translated equal-size intervals [0 + xir/10, tt/2 + xtv/10], 6a and [— 7r/4 + 



, 7r/4 + x7r/20], 6b with the different colors, going from blue-green to red, showing x 

3,4,5. 




FIG. 7. AS v n f° r intervals of various sizes [tt/4 — xtt/20, tt/2 + xtt/20], with the different colors, 
going from red to blue-green, showing x = 1, 2, 3, 4, 5. 



VI.2. The entanglement from the deformation operator 



Here we examine the second Renyi entropy that results from quenching with the full 
deformation operator. The expression for S2 in (5.3) may be written as 



S_ ovac 1 chare riins 
2 — 2 + 2 -+- 2 ) 



(6.11) 



where S 2 ac is the Renyi entropy of the vacuum in (5.4), S 2 are is the additional Renyi en- 
tropy added by the bare twist in (5A_), and 5™ s is the new contribution from the super- 
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charge and spin field in the covering space. We write the three pieces as 



tjivac 



gbare 



s ms 



log 
log 

log 



R 2 



'23 



Z2Z3 
16(77] 



1 + Vv\' 



1 — 774 


10 + |1 + 773 


10 


2io 


n 


5 
4 



(6.12a) 
(6.12b) 

(6.12c) 



Using (|6.4|), we can rewrite the last two terms as 
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(6.13a) 
(6.13b) 



and AS 2 = S% are + S^ s . We plot this for the interval [n/2, 3tt/4] in Figure |4b| and, as men- 
tioned, plots analogous to those in Figures J6aj|6bj and|7]look very similar. Let us note that 
S^ s is distinguished from the entanglement of the vacuum and of the bare twist in that 
the left and right contributions do not directly factorize. That being said, one finds that 
gins g^rj en j y S a superposition principle with respect to the left and right contributions 
as long as left- and right-moving excitations do not simultaneously contribute to the entropy. 
To illustrate consider in the limit of vanishing e. In particular, the peak value of 
when left or right excitations are separately in the interval is given by 



limSr 



log 2, 



(6.14) 



sing. cont. 

and is finite. When both left and right excitations contribute, however, Sf* is given by 



limSf 3 



4 log 2, 



(6.15) 



both cont. 



which serves as a finite upper bound on S^ s ; however, it is not simply twice the con- 
tribution from left- or right-movers separately. Aside from the failure of left and right 
contributions to factorize, is further distinguished in its finiteness. In contrast, S^ 316 
represents a huge (actually divergent) amount of entanglement. This suggests that it is 
the twisting part of the deformation operator that predominantly contributes to thermal- 
ization. This is consistent with the point of view advocated in Il3"6l438ll , which focuses on 
the effect of the twist operator. 



VII. CONCLUSION 

We have seen how to analytically compute entanglement entropies for arbitrary spatial 
intervals of the D1D5 CFT following a local quench of the system by an exactly marginal 
deformation operator, which contains the 2-twist operator and moves the theory toward 
the supergravity regime in its moduli space. For the insertions of the bare twist operators 
we were able to compute all of the Renyi entropies S n and the von Neumann entropy 



25 



while for the full deformation operators we were only able to compute the S n for n > 
2. We could understand the qualitative behavior of the entropies as arising from pairs 
of entangled particles generated at the quench event and propagating from there at the 
speed of light. 

This process does not lead to thermalization of the system, as we saw in the discus- 



sion below (6.9), but from the point of view of the interval these particles should appear 
as thermal radiation in accord with their nonzero entropy The nonzero entropy does not 
itself, of course, guarantee the radiation is thermal, but we expect it is for several reasons. 
First we note that in general global quenches, which can also be understood as producing 
pairs of entangled particles ||53|, do lead to thermalized systems ||T0llTTll50| . In our case of 
a local quench, the positive entanglement entropy is apparently due to the entanglement 
of pairs of localized excitations, one member of which is inside the interval and one out- 
side, and so is traced over. This is a situation familiar from Hawking radiation: it appears 
thermal to observers outside the horizon. 

A toy model illustrates this in more detail. Consider two harmonic oscillators with 
operators {a, a)} and {b, tf} in the entangled state \ip) oc e Aatfot |0). Such a squeezed state is 
the general result of a Bogolyubov transformation on the operator algebras and appears 



in our case as (3.1) above. If we then trace over, say, the b oscillator then the resulting 



reduced density matrix is given by 

p a oc ^ \\\ 2n \n) a (n\ a . (7.1) 

n 

Now we can compare to a thermal density matrix for the a oscillator using ( |1.1| >: 

p th ex e"^ cx 5>-*"» a H a , (7-2) 

n 

and we see that p a is a thermal density matrix at inverse temperature = — - log |A| 2 . It is 
straightforward to compute the entropies of the state p a as a function of A (or equivalently 
(3) but we do not need the explicit formulas here. 

This simple calculation illustrates how thermal density matrices, with their associated 
entropies, arise from the process we are considering. Of course, the trace we perform in 
the CFT, over the exterior of a spatial interval, is not directly analogous to this simple 
case. It seems to be difficult to perform such a spatial trace directly on the state after the 
quench, given by ( |3.1| >, which is why we pursued a technique here that makes extensive 
use of the powerful conf ormal symmetry of the system. It would be interesting to pursue 
that direct approach and also to trace over different classes of subsystems, corresponding 
to different coarse-grainings. 

Ultimately we would like to understand thermalization in this system, with general 
subsystems characterized by reduced density matrices of the kind in ( 1.1) . Here we have 



only studied an individual event involving a small sector of the full theory, which is a 
system of N]_ Dl and iV 5 D5 branes with N\ and N 5 potentially large. This has been in the 
spirit of time-dependent perturbation theory of a weakly interacting system: individual 
interactions are treated separately and the cumulative effect of many independent inter- 
actions is put together at the end. We can imagine how the thermal radiation produced 
by many independent local quenches can ultimately lead to a thermalized system, but a 
full, careful treatment remains to be done. 
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We would like to compare the entropies produced from such a process to those of the 
system in equilibrium at some finite temperature, not just the high-temperature limit we 
considered above. Computing the finite-temperature entropies would involve comput- 
ing two-point functions of twist operators in a domain with both space and time period- 
ically identified, i.e., a torus. This is somewhat challenging as the answer would depend 
on the full operator content [62J and one cannot uniformize to the plane. It has been 
carried out for a single fermion |67| and we hope to address this for the D1D5 system in 
future work. It should also be possible to learn more information about the state after 
the quench by judiciously studying the whole set of Renyi entropies that can be obtained 
from (4.20) and ( 4.23| > (or ( |5.1| > for the bare twists), rather than just n = 1 or 2 as we did 
here, since they collectively contain significantly more information. Our results (partic- 
ularly for the bare, nonsupersymmetric twist operators) may be relevant to and possibly 
subject to verification by the local quenches studied in condensed matter systems, e.g., 
|5Tll68]470| . However, there are a number of issues in making such a comparison, since 
our CFT and local quench are both highly specific, and we have not seriously attempted 
to do so. 

Finally a few words on the bulk description of the process we studied. The Dl-branes 
are wrapped in the D5-branes, which are localized in the transverse asymptotically 4+1 
dimensional Minkowski space. Initially the branes are in a stationary state corresponding 
to the ground state of the D1D5 orbifold CFT. We then imagine a sudden interaction with 
an external field in the transverse space, which weakly deforms the D1D5 system by the 
exactly marginal twist operator (2.4) and supplies the necessary energy for thermalization 
to begin. This is admittedly just a sketch that we hope to improve. 
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